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Abstract 

We investigate the phase structure of lattice QCD for general number of flavors 
Np. Based on numerical data combined with the results of the perturbation theory 
we propose the following picture: When Np > 17, there is only one IR fixed point at 
vanishing gauge coupling, i.e., the theory in the continuum limit is trivial. On the other 
hand, when 16 > Np > 7, there is a non-trivial fixed point. Therefore, the theory is 
non-trivial with anomalous dimensions, however, without quark confinement. Theories 
which satisfy both quark confinement and spontaneous chiral symmetry breaking in 
the continuum limit exist only for Np < 6. 



*' Talk presented by K. Kanaya at the 1997 Yukawa International Seminar (YKIS'97) on "Non- 
Perturbative QCD — Structure of the QCD Vacuum — ", YITP, Kyoto, Japan, 2-12 Dec. 1997. To be 
published in the proceedings [Prog. Theor. Phys. Suppl.]. 



1 typeset using pTpTfXsty <ver.0.8> 



1. Introduction 



At low temperatures, QCD with a small number of light quarks is known to have three 
characteristic properties of asymptotic freedom, quark confinement and spontaneous break- 
down of chiral symmetry. Among them, asymptotic freedom is reliably investigated by the 
perturbation theory. It is well-known that when the number of flavors Np exceeds 16 1, 
the asymptotic freedom is lost. Then, a question which naturally arises is whether there 
is a constraint on Np for quark confinement and/or the spontaneous breakdown of chiral 
symmetry. Because these properties are essentially non-perturbative in nature, we have to 
carry out a non-perturbative investigation. 

In this paper, we study the Np dependence of the QCD vacuum by numerical simulations 
on the lattice. Because the Wilson fermion formalism is the only formalism for lattice quarks 
that preserves manifest flavor symmetry, we use the Wilson fermion formalism in our study 
of Np dependence in QCD. Based on the numerical data obtained for Np from 2 up to 
300, we elucidate the phase structure of QCD as a function of Np. We conjecture the 
following picture: When Np > 17, there is only one IR fixed point at vanishing gauge 
coupling. Therefore, the theory in the continuum limit is trivial. On the other hand, when 
16 > Np > 7, there is a non-trivial IR fixed point. Therefore, for this case, the theory 
is non-trivial with anomalous dimensions, however, without quark confinement. Theories 
which satisfy both quark confinement and spontaneous chiral symmetry breaking in the 
continuum limit exist only for Np < 6. Our results are in part reported in Refs. |l|) and |^). 

This paper is organized as follows. Before going into details of our results obtained from 
numerical simulations, in Sec. ^|, we first summarize results and speculations from a pertur- 
bative study of the renormalization group beta function, and then present our conjecture 
for the Np dependence of the beta function, based on our numerical data obtained on the 
lattice. In Sec. || our lattice model and simulation parameters are described. In order to in- 
vestigate the Np dependence of the QCD vacuum in the continuum limit, we study the phase 
structure of lattice QCD at zero temperature for general number of flavors, in particular, for 
Np > 7. When the phase diagram becomes clear, we are able to see the nature of the QCD 
vacuum in the continuum limit, and eventually answer the question about the condition on 
Np for confinement and spontaneous breakdown of the chiral symmetry. Results for the 
phase structure in the strong coupling limit are given in Sec. |j. Sec. [5] is devoted to the 
phase structure off the strong coupling limit. We then discuss the nature of the deconfining 
phase in Sec. ||. Our conclusions and conjectures concerning the phase structure of QCD for 
general number of flavors are summarized in Sec. |7|. 
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b) 

Fig. 1. Renormalization group beta function, (a) Conjecture by Banks and ZaksEP assuming con- 
finement in the strong coupling limit for all Np. (b) Our conjecture deduced from the results 
of lattice simulations. 



§2. Beta function 

The perturbative beta function of QCD is universal up to two-loop order: (3(g) = —b g 3 — 
b x g 5 + ■■■ .where b = ^ (ll - |JV>) and b x = (l02 - fNp) . The coefficient b 

changes its sign at Np = 16|. Therefore, for Np > 17, asymptotic freedom is lost and the 
IR limit is free, i.e. we have no confinement. It is well-known that the second coefficient 
bi changes its sign already at Np ~ 8.05. The values of the two coefficients imply that for 
N F ~ 9, the IR fixed point appears at g ~ 3. However, this value of the gauge coupling 
constant is too large to trust the result from the perturbation theory. In contrast to this 
case, for Np ~ 16, the IR fixed point appears in a perturbative region. Therefore, it is 
plausible that the beta function has a non-trivial IR fixed point for N' < Np < 16 with 
some N' < 16. When such an IR fixed point exists, the coupling constant cannot become 
arbitrarily large in the IR region — this implies that quarks are not confined. 

A pioneering study on the Np dependence of the QCD vacuum was made by Banks and 
Zaks in 19820*. Based on an early result obtained on the lattice they assumed that the 
beta function is negative in the strong coupling limit for all Np. Using the perturbative 
results for N' < Np < 16 mentioned above, they conjectured Fig. 1(a) as the simplest Np 
dependence of the beta function, and studied the phase structure of QCD based on this Np 
dependence of the beta function.0 Because of an additional non-trivial UV fixed point for 

*' For different approaches, see Refs. @),§) and g). When the quarks are confined, the sigma model might 
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Np > N', their conjecture for the phase structure is complicated. 

The assumption of confinement in the strong coupling limit plays an essential role in the 
arguments of Banks and Zaks. Here, a quite important fact to be noted is that the lattice 
study cited by them is a study of pure gauge theory without quarks &. 

When dynamical quarks are coupled, the vacuum structure become much complicated. 
In this case, we can formally rewrite the theory in terms of mesonic and baryonic fields [see, 
for example, Ref. |9])]. When mesons and baryons are physically relevant dynamical variables, 
i.e. when quarks are confined, the effective action for hadrons is useful to study the phase 
structure of the system. Oppositely, when the resulting effective action for hadrons is a 
well-defined action of weakly interacting particles, it is plausible that quarks are confined. 
In the strong coupling limit, computation of the effective theory for hadrons are carried out 
either in the large N c limit i\ using a meanfield approximation (1/d expansion) or 
using a heavy quark mass expansionE^'H*. In these cases, the effective action for mesons 
and baryons leads to a world with spontaneously broken chiral symmetry. Here, in contrast 
to these cases, we are interested in the case where N c = 3 and quarks are light. In this case, 
the effective action is quite complicated for general Np; diagrams with quark loops, which 
are suppressed at large N c , large d or large m q , become important. Therefore, in the realistic 
cases, it is difficult to analytically deduce the vacuum structure as a function of Np even in 
the strong coupling limit. 

Therefore, in our previous work^, we performed a series of numerical simulations of QCD 
in the strong coupling limit for various Np, using the Wilson fermion formalism for lattice 
quarks. We found that, when Np > 7, light quarks are deconfined and chiral symmetry is 
restored at zero temperature even in the strong coupling limit. 

We extend this study to weaker couplings and to larger Np. Based on the numerical 
results obtained on the lattice combined with the results of the perturbation theory (see the 
following sections for details), we conjecture Fig. 1(b) for the Np dependence of the beta 
function: When Np is smaller than a critical value, the beta function is negative for all 
values of g. Quarks are confined and the chiral symmetry is spontaneously broken at zero 
temperature. We find that the critical number of flavors is 6. (Corresponding critical Np 
is 2 for the case N c = 20>. ) On the other hand, when N F is equal or larger than 17, we 
conjecture that the beta function is positive for all g, in contrast to the conjecture by Banks 
and Zaks shown in Fig. 1(a). The theory is trivial in this case. When Np is between 7 and 
16, the beta function changes sign from negative to positive with increasing g. Therefore, 
the theory has a non-trivial IR fixed point. 



be a reasonable effective theory of QCD. The ./Vp-dependence of the vacuum in a sigma model is studied in 
Ref. 0). 



4 



§3. Model and simulation parameters 
Lattice QCD is denned on a hypercubic lattice in 4-dimensional Euclidian space with 



lattice spacing a. Formulation of lattice QCD is described in Ref. [14] ) . 

In this study, we use the standard one plaquette gauge action combined with the standard 
Wilson quark action: S = S gauge + S quar k with 

Sgauge = — ^ Tr {^Ux^Ux+fi^U^^p^U^p^ , (3'1) 



Squark = £ $M ~ K j ( 1 - 7 M ) U x ^ f x+(k + tfl+pfr+lpWl^l 



(3-2) 



where = 6/g 2 with (7 the bare gauge coupling constant.^ The index / (= l,---,Np) 
is for the flavors. We assume that all Np quarks are degenerate with the bare mass m , 
which is related to the hopping parameter K by itlq = (l/2a)(l/K — 1/K C ), where a is the 
lattice spacing and K c — 1/8 is the point where the bare mass vanishes. The point K = 
corresponds to infinite quark mass. 

In the formalism of Wilson quarks on the lattice, the flavor symmetry as well as C, 
P and T symmetries are exactly satisfied on the lattice with a finite a. However, chiral 
symmetry is explicitly broken even at m = 0. [For details about the effects of broken 



chiral symmetry, see Ref. [16]) and references cited there.] However, through a perturbative 
study of Ward identities using the Wilson fermion action0, it is shown that the effects of 
the broken chiral symmetry can be absorbed by appropriate renormalizations, including an 
additive renormalization of the quark mass, i.e. by a shift of K c as a function of (3. 
We define the current quark mass in terms of an axial vector Ward identity: 

2m q (0\P\7r(p = 0)) = -m n ( | A A | 7r(p = 0) ) (3-3) 

where P is the pseudoscalar density, A4 the fourth component of the local axial vector 
current, and m n the pion (screening) massli^r A multiplicative renormalization factor 
for the axial current, which is not important in this study, is absorbed into the definition 
of the quark mass. Numerical simulations show that the value of m q does not depend on 
whether the system is in the high or the low temperature phase when (3 is large; (3 > 5.5 for 
N F = 2©. 

We define the chiral limit K c by the condition m q = 0. Numerical simulations show that, 
for small A?^, K c (0) forms a smooth curve connecting 1/8 at (3 = 00 and 1/4 at (3 — 0. When 
we further define the pion decay constant f n by ( | A4 1 7r(p = 0) ) = m^f^, either m n = 



*) 



Many flavor QCD using the staggered fermion formalism is studied in Ref. [La) 
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or fn = is satisfied in the chiral limit. The chiral symmetry is restored when f w = in the 
chiral limit, while = if the chiral symmetry is spontaneously broken. 

For QCD with a small Np, the deconfinement transition is expected at high temperatures 
0). Therefore, we have to study the temperature dependence of the results to distinguish 
between the deconfinement transition due to high temperatures and a bulk transition due 
to the effects of many flavors. On a lattice with N t sites in the Euclidian time direction, the 
temperature is given by T = 1/N t a. When the beta function is negative, as in the case of 
small Np, the lattice spacing a is a decreasing function of (3 = 6/g 2 . In this case, T increases 
with increasing (3 for a fixed finite N t . On the other hand, when the beta function is positive, 
T decreases with increasing (3. Therefore, in order to distinguish between the deconfinement 
transition due to high temperatures and a bulk transition for general Np, we have to vary 
N t . 

We perform simulations on lattices 8 2 x 10 x N t (N t = 4, 6, and 8), 16 2 x 24 x N t (N t = 16) 
and 18 2 x 24 x N t (N t = 18). We vary N F from 2 to 300. For each Np, we study the phase 
structure in the coupling parameter space (f3,K). We adopt an anti-periodic boundary 
condition for quarks in the t direction and periodic boundary conditions otherwise. We use 
the hybrid R algorithm© for the generation of gauge configurations.0 Statistical errors are 
estimated by the jack-knife method. 

It should be noted that, in QCD with dynamical quarks, there are no order parameters 
for quark confinement. We discuss about confinement by comparing the screening pion mass 
and the lowest Matsubara frequency, and, simultaneously, consulting the values of plaquette 
and the Polyakov loop. In the followings, we call the pion screening mass simply the pion 
mass, and similarly for the quark mass. 



§4. Strong coupling limit 

We first study the strong coupling limit, (3 = 0, where we expect the strongest confining 
force. In a previous work©, we have shown that, when Np > 7, light quarks are deconfined 
and chiral symmetry is restored at zero temperature even in the strong coupling limit: For 
Np < 6, we have only one confining phase from the heavy quark limit K = up to the chiral 
limit K c = 0.25. On the other hand, for Np > 7, we find a strong first order transition at 
K = K d < K c . When quarks are heavy (K < K d ), both plaquette and the Polyakov loop 
are small, and m n satisfies the PCAC relation m 2 oc m q . Therefore, quarks are confined and 

*> The R algorithm has discretization errors of 0(NpAr 2 ) for the step size At of a molecular dynamic 
evolution. As Np increases we have to decrease At, such as At =0.0025 for Np = 240, to reduce the errors. 
We check that the errors are sufficiently small with our choices of At for typical cases. 
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Fig. 2. Plaquette at (5 = as a function of 1/K. (a) For various N F at N t = 4. (b) For N F = 240 
at N t = 4, 8 and 16. 1/K d ~ 8.1(1) for N t = 4 and 1/K d = 7.8(2) for N t = 8 and 16. 
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Fig. 3. The transition point 1/Kd at (3 = versus Np for Nt = A and iVt > 8. For clarity, data at 
Nt = 8 for iVp = 300 is slightly shifted to a larger Np in the figure. 



the chiral symmetry is spontaneously broken in this phase. We find that m q in the confining 
phase is non-zero at the transition point K d , i.e., the chiral limit does not belong to this 
phase. On the other hand, when quarks are light (K > K d ), plaquette and the Polyakov 
loop are large. In this phase, m v remains large in the chiral limit and is almost equal to twice 
the lowest Matsubara frequency ir/N t . This implies that the pion state is a free two-quark 
state and, therefore, quarks are not confined in this phase. The pion mass is nearly equal to 
the scalar meson mass, and the rho meson mass to the axial vector meson mass. The chiral 
symmetry is also manifest within corrections due to finite lattice spacing. 

We now extend the study to larger Np- In Fig. 2(a), we show the results of plaquette 
at Nt = 4 for Np = 7 - 300. Clear first order transition can be seen at K below 0.25. We 
then study the N t dependence of the results, as shown in Fig. 2(b) for the case of Np = 240. 



Results for Np = 7 is given in Ref. ^3|) . We find that, although the transition point shows 
a slight shift to smaller 1/K when we decrease T = 1/N t a from Nt = 4 to 8, it stays at the 
same point for N t > 8. The chiral limit, 1/K = 4, cannot be achieved in the confining phase 
even at T = 0. Therefore, we conclude that the transition is a bulk transition (transition 
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Fig. 4. The phase structure for Np < 6; (a) at finite temperatures, and (b) at zero temperature. 
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Fig. 5. Phase diagram for (a) Np = 12 and (b) Np = 18. Dark shaded lines represent our 
conjecture for the bulk transition line in the limit N t = oo. 



that exists at zero-temperature). From these results, we obtain the phase diagram at (5 = 
shown in Fig. 3. 



§5. Phase structure at (3 > 

Let us now study the case (3 > 0. On a lattice with a fixed finite Nt, we have the finite 
temperature deconfining transition at finite (3 when the quark mass is sufficiently heavy 
(K pa 0). For iVp < 6, our previous study shows that the finite temperature transition line 
crosses the K c line at finite /3@>. A schematic diagram of the phase structure for this case 
is shown in Fig. 4.0 The location of the finite temperature transition line moves towards 
a larger (3 as N t is increased. In the limit N t = oo, the finite temperature transition line 
will shift to (3 = oo so that only the confining phase is realized for iYV < 6 at T = 0. See 
Fig. 4(b). There, we have a monotonous RG flow towards the strong coupling limit on the 
K c line. 

When Np > 7, we have seen that there exist a bulk deconfining transition in the strong 
coupling limit. We have to clarify the relation between the finite temperature transition at 
small K and the bulk transition at /3 — 0. 

*) See Ref. |l^) and references cited there for discussions about complexities due to the lack of chiral 
symmetry with Wilson fermions. 
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Fig. 6. Phase diagram for Np = 240. 

At (3 — 0, the location of the bulk deconfining transition Kd moves toward larger K with 
decreasing N F . From Fig. 3, we find 1/K d ~ 8.5, 8, 7.6, 7.2, 6.1, and 5.0 for Np = 300, 240, 
180, 120, 60, and 18, respectively, at /3 — 0. Recall that 1/K — 8 is the point for massless 
free quarks. 

We perform simulations at f3 « 0.0 - 6.0 varying Np from 2 to 300. From the numerical 
data, we obtain the phase diagrams sown in Fig. 5 for Np = 12 and 18, and in Fig. |5] for Np = 
240. When Nt — 4 or 8, the phase boundary line between confining and deconfining phase 
bends down at finite (3 due to the finite temperature phase transition of the confining phase. 
The dark shaded lines are our conjecture for the phase boundary between the confining phase 
and the deconfining phase at zero temperature (N t = oo). 

In Figs. 5 and |6|, the dashed lines in the deconfining phase represents the chiral limit, 
m q = 0. This line also corresponds to the minimum point of m^. We have checked that, 
for the case Np = 240, the quark propagator in the Landau gauge actually shows the chiral 
symmetry, 75^(^)75 = —G(z), on the m q = line. For A^p>240, we find that the m q = 
line in the deconfining phase, which starts at 1/K = 8 in the weak coupling limit (3 = 00, 
reaches the strong coupling limit, as shown in Fig. |6] for Np = 240. For a smaller iVV < 100, 
because the bulk transition line in the strong coupling region shifts toward larger K with 
decreasing Np, the m q = line in the deconfining phase hits the line at finite (3. For 
example, in the case Np = 18, it hits at (3 = 4.0 - 4.5, as shown in Fig. 5(b). 



§6. Nature of the deconfining phase for large Np 

In order to study the nature of the deconfining phase, we first study the case of large iVV, 
because the range of the deconfining phase is wider at larger Np (see Fig. 3) and therefore 
the parameter dependence of physical quantities becomes easier to study. In particular, for 
Np > 240, we are able to study the chiral limit in the deconfining phase at any values of (3 
(see Fig. ^j. Because of this, we first intensively investigate the cases N F = 240 and 300, 
and then decrease N F . 
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Fig. 7. Results of ml and 2m q versus 1/K at N t = 4: (a) N F = 240 at various (5. (b) iV F = 18 - 
300 at f3 = 6.0 and 0.0. 



6.1. N F = 240 and 300 

Fig. 7(a) shows the results of ml and 2m q in the deconfining phase for Np = 240 at 
(3 = 0.0, 2.0, 4.5, 6.0, and 100.0 on the N t = 4 lattice. A very striking fact is that the shape 
of ml and 2m q as a function of 1/K only slightly changes for 1/K < 8 when the value of f3 
decreases from oo down 0. [Note that (3 = oo corresponds to the case of free quarks. Due to 
the doublers, m q defined by (|3-3|) shows a non-trivial 1/K dependence at 1/K < 8 also for a 
free quark.] Only the position of the local minimum of at 1/K ~ 8, which corresponds 
to the vanishing point of m q , slightly shifts toward smaller 1/K. We obtain similar results 
also for N t = 8. This suggests that quarks are almost free down to (3 = in the deconfining 
phase. 

The results for Np = 300 are essentially the same as those for Np = 240, except for very 
small shifts of the transition point and the minimum point of m\. 

6.2. Renormalization group flow for Np = 240 

From the perturbation theory, the m q = point at (3 = oo is a trivial IR fixed point for 
Np > 17. The phase diagram shown in Fig. ^ suggests that there are no other fixed points 
on the m q = line at finite (3. In order to confirm this, we investigate the direction of the 
renormalization group (RG) flow along the m q = line for Np = 240, using a Monte Carlo 
renormalization group (MCRG) method. 

We make a block transformation for a change of scale factor 2, and estimate the quantity 
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Aft = ft {2a) — /3(a): We generate configurations on an 8 4 lattice on the m q = points at 
(3 = and 6.0 and make twice blockings. We also generate configurations on a 4 4 lattice and 
make once a blocking. Then we calculate Aft by matching the value of the plaquette at each 
step.0 From a matching using our data, we obtain Af3 ~ 6.5 at ft = and 10.5 at ft = 6.0. 
The value obtained from the two-loop perturbation theory is Af3 ~ 8.8 at ft = 6.0. The 
signs are the same and the magnitudes are comparable. This suggests that the direction of 
the RG flow on the m q = line at (3 = and 6.0 is the same as that at ft = oo. This further 
suggests that there are no fixed points at finite ft. All of the above imply that the theory is 
trivial for N F = 240. 

6.3. 240 > N F > 17 

Now we decrease Np from 240. As discussed in Sec. [5], the area of the deconfining phase 
decreases with decreasing Np in the strong coupling region. However, except for this shift of 
the bulk transition point, the shape and the values of and m q are quite similar when we 
vary Np from 300 down to 17. Fig. 7(b) shows m% and 2m q for Np = 18 - 300 at (3 = 6.0 
and 0. At (3 = 6.0, the shapes of are almost identical to each other, except for a small 
shift toward smaller 1/K as Np is decreased. 

These facts suggest that, for Np > 17, the nature of physical quantities in the deconfining 
phase is almost identical to that observed for Np = 240 and 300, i.e., quarks are almost free. 
Therefore, we conjecture that the direction of the RG flow at (3 < 6 is identical to the case 
Np = 240. Combining this with the perturbative result that, for Np > 17, (3 = oo is the IR 
fixed point, we conjecture that the RG flow along the massless quark line in the deconfining 
phase is uniformly directing to smaller ft, as in the case of N F = 240. When this is the case, 
the theory has only one IR fixed point at ft = oo — the theory is trivial for Np > 17. 

6.4. 16 > Np > 7 

As shown in Sec. quark confinement is lost for Np > 7 at ft = 0. We have intensively 
simulated the cases Np = 12 and 7. The phase diagram for Np = 12 is shown in Fig. 5(a). 
At ft < 6.0, the gross feature of the phase diagrams for Np = 7 and 12 is quite similar to the 
case Np = 18 shown in Fig. 5(b). Physical quantities at ft <6.0 also show behavior similar 
to that shown in Fig. 7 for Np > 17. Therefore, we consider it probable that the direction 
of the RG flow in the deconfining phase at small ft is towards a larger ft as in the case of 

*' It is known for the pure SU(3) gauge theory, that one has to make a more careful analysis using several 
types of Wilson loop with many blocking steps to extract a precise value of A(3. We reserve elaboration of 
this point and a fine tuning of 1/K at each (3 for future works. For Np = 240, because the velocity of the 
RG flow is large, we will be able to obtain the sign and an approximate value of A(3 by a simple matching. 
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Fig. 8. Phase structure and the RG flow at T = for general number of flavors. Thin dashed lines 
in the phase diagrams for Np > 7 represent the location of the chiral limit for Np < 6 as a 
guide for eyes. 

Np > 17. 

On the other hand, when Np < 16, the theory is asymptotic free. Therefore, the RG 
flow at (3 = oo is opposite to that for Np > 17. This means that we have an IR fixed point 
somewhere at a finite value of (3. The continuum limit is governed by this IR fixed point. 

On the m q = line around /3 = 4.5, we find that m n is roughly twice the lowest Matsubara 
frequency. This implies that the quarks are not confined and almost free. Therefore, the 
anomalous dimensions are small at the IR fixed point, suggesting that the IR fixed point 
locates at finite (3. In order to pin down the position of the IR fixed point, we have to make 
a detailed study of RG flows in a wider coupling parameter space. We reserve these studies 
for future works. 

In summary, we conjecture that, for 16 > Np > 7, we have an IR fixed point at finite (3. 
The theory in the continuum limit is a non-trivial theory with small anomalous dimensions, 
however, without confinement. 



§7. Conclusions 

Based on our numerical data combined with the results of the perturbation theory, we 
propose Fig. 8 for the phase structure and the RG flow at T = for general number of 
flavors. 
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When Np < 6, quarks are confined for any values of the current quark mass and (3. 
The chiral limit K c (/3), where the current quark mass m q vanishes, belongs to the confining 
phase. 

When Np > 7, there is no chiral limit in the confining phase. There is a line of a first 
order phase transition from the confining phase to a deconfining phase at a finite current 
quark mass for all values of f3. The chiral limit exists only in the deconfining phase. 

In order to study the nature of the deconfining phase for large Np, we first increase Np, 
because the region of the deconfining phase becomes larger with larger Np. For Np = 240, a 
MCRG study along the massless quark line shows that the direction and the magnitude of the 
RG flow at small f3 is consistent with the result of the perturbation theory, suggesting that 
the direction of the RG flow is identical for all values of (3, as shown in Fig. 8. Therefore, 
there is only one IR fixed point at (3 = oo, so that the theory in the continuum limit is 
trivial. We also find that the independences of pion screening mass and current quark 
mass in the deconfining phase are almost identical to those of free Wilson quarks. Also for 
smaller values of Np > 7, the general features of the phase diagram and the i^-dependence of 
physical quantities are quite similar to the case of large Np ~ 240. Therefore, we conjecture 
that the direction of the RG flow is the same as in the case of N F = 240, at least for f3 < 6.0 
we simulate. 

For Np > 17, because the RG flow in the weak coupling limit has the same direction due 
to the lack of asymptotic freedom, we conjecture a uniform RG flow towards the trivial IR 
fixed point, as in the case of Np > 240. Thus the theory is trivial for Np > 17. 

When 7 < Np < 16, because the theory is asymptotically free, we expect a non-trivial 
IR fixed point at finite f3. In this case, the theory in the continuum limit is a non-trivial 
theory with anomalous dimensions, however, without confinement. 
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